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This paper presents some applications of the vortex-lattice theory to the 
preliminary aerodynamic design and analysis of subsonic aircraft. These methods 
include the Rockwell-Tulinius vortex-lattice theory for estimating aerodynamic 

I characteristics, a Trefftz plane optimization procedure for determining the span 
loads for minimum induced drag, and a modification of the TrefTtz plane proce- 
dure to estimate the Induced drag for specified span loads. The fi ;t two 
methods are used to aerodynamically design aircraft planforms, twists, and cam- 
bers, and the latter method is used to estimate the drag .'.'or components such as 
flaps and control surfaces. 

Results from the theories for predicting lift and pitching moment, drag dur 
to lift, and the drag of control surfaces are compared with experimental data. 
The data were obtained on a general aviation model with flaps and a close- 

' coupled canard-wing model. 
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In the preliminary stages of aircraft design, it is necessary that the 
designer have valid estimates of aircraft aerodynamics, particularly lift, drag, 
and pitching moments. Lift and pitching moment are required to size the plan- 
forms (wing, tall, and canard) and locate them with respect to a moment center, 
usually a desired aircraft center of gravity, for trimmed lift requirements and 
stability margins. Skin friction, form, and Induced drags must be estimated 
and minimized for best performance. Many theoretical methods involving various 
levels of complexity have been developed which estimate these characteristics to 
varying degrees of accuracy. The preliminary designer, however, wants methods 
that are fast, reasonably accurate, and easy to use so that changes in aircraft 
configuration can be easily assessed. Once the overall configuration geometry 
is defined, he may wish to use some of the more highly sophisticated methods to 
refine his estimates before beginning experimental verification of the design. 
This paper will address applications of easy-to-use methods appropriate at the 
preliminary design stage; these methods Include the Rockwell-Tulinius vortex- 
lattice theory for estimating aerodynamic characteristics, a Trefftz plane opti- 
mization procedure for determining the span loads for minimum induced drag, and 
a modification of the Trefftz plane procedure to estimate the induced drag for 
specified span loads. 
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SYMBOLS 

aspect ratio, b^/S 
span 

drag coefficient 
induced drag coefficient 
minimum drag coefficient 
lift coefficient 
lift-curve slope 

pitching-moment coefficient 
normal-force coefficient 
chord 

average chord 
mean aerodynamic chord 
section li't coefficient 
section normal-force coefficient 

y-component of influence function for pair of trailing vortex legs 

z-component of influence function for pair of trailing vortex legs 

induced-drag efficiency parameter, i^^ 

vertical separation -between canard and wing 
y-component of normal unit vector 

z-component of normal unit vector 

dynamic pressure 
wing area 

incremental section width (from ref. 1) 
tangent unit vector spanwise component 
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tangent unit vector vertical 


free-stream velocity 

W 

downwash velocity 

X,Y,Z 

axis system 

x,y 

distance along X- and Y-axes 

X 

moment-center location 

me 


a 

angle of attack 

r 

section circulation 


flap deflection 

n 

fraction of semispan, y/b/2 

p 

density 

Subscripts: 

cp 

center of pressure 

d 

design 


indices 

max 

maximum 


DISCUSSION 


Prediction of Lift and Pitching Moment 

The Rockweil-Tulinius unified vortex-lattice theory (refs. 1 and 2) can be 
used to predict static and rotary stability derivatives for configurations with 
multiple lifting surfaces of arbitrary shape. It can also compute the section 
and total configuration forces and moments for arbitrary planform geometries 
with twist and camber. This method, as programed, is fast, easy to use, and 
fairly accurate. 

The agreement between this theory and experimental data for the lift of a 
simplified general aviation model is shown in figure 1. The model has a 
straight untapered wing using the NASA GA(W)-1 airfoil section (refs. 3 and 4) 
and had 2° of twist (washout) from the root to the tip. The model body was a 
flat-sided ellipse. For the theoretical calculations, the fuselage was modeled 
as a flat plate and the wing as a camber line with twist. Agreement between the 
estimated and the experimental C was quite good at low angles of attack 

prior to flow separation which occurred at a = 4°. 
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Also shown in figure 1 are the theoretical drag polars for 0-percent and 
100-percent leading-edge suction as given by the equations 
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for 100-percent leading-edge suction and 



for 0-percent leading-edge suction. 


The value for 

D,o 


was obtained from the 


experimental data. These curves for 100-percent and 0-percent leading-edge 
suction represent the best and worst possible drag polars, respectively, for a 
given configuration. The leading-edge radius and/or camber design should pro- 
duce data that are as close to the 100-percent suction polar as possible. Near- 
field analyses are required to minimize viscous and separated flow effects to 
approach the 100-percent suction polar. For this case, the data show that the 
design was close to the 100-percent suction polar up to = 1.2. 


Two-dimensional separation can be delayed and minimized for moderate angles 
of attack by proper planform shaping, camber design, and leading-edge-radius 
selection. However, at large angles of attack, the viscous form drag must be 
reduced by taking advantage of interfering flow fields of adjacent surfaces, 
vortex flows, or induced propulsion effects. Examples of applications of this 
appiuach for reducing viscous form drag due to lift are shown later in this 
paper. 


Figure 2 is a sketch of a close-coupled canard model tested in the Langley 
V/STOL tunnel to investigate the effects of propulsion on stability at high 
angles of attack. A similar unpowered model was tested in the Langley high- 
speed 7- by 10-foot tunnel by Blair B. Gloss (ref. 5) to determine the effect of 

vortex lift on performance, especially C . The wings and canards of both 

,max 

models had symmetrical circular arc airfoil sections. Also, strakes were uti- 
lized in both tests to produce vortex lift at the higher angles of attack. The 
agreement between theory and data of Gloss (fig. 3) is good over the linear 
range of the data for the wing and the wing-canard configurations. The method 
does not predict the additional vortex lift and resulting pitching moment when 
the strake is present. 

This method was used to establish a moment center for a wing-canard model 
to give a stability margin at low of -5 percent ” 0.05^ prior to 

testing in the V/STOL tunnel. The data, shown in figure 4, indicate a value of 

9C /9C^ of about 0.06 to 0.05 at low C. , which agrees well with the predi ted 
m/ L L 

value. 
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Prediction of Minimum Induced Drag 

Once the plai.Jorms of a configuration have been sized and located to meet 
lift and stability requirements, it is necessary to compute che optimum span 
loads for minimum induced drag for the interfering planforms. The expression 
for the induced drag was developed by using an equivalent lifting-line Trefftz 
plane approach of reference 1 and is illustrated by the following sketch and 
equation: 
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By utilizing the method of Lagrangian multipliers with the induced-drag equation, 
the span loads for minimum C . may be calculated while constraining C. and 

U y X L< 

C to desired values. This procedure was programed by Tulinius and Gloss, and 
m 

the results are given in reference 2. The input for this program consists of 

the basic planform geometry, as in the Rockwell-Tullnlus method, along with the 

desired (x/c) distribution. The (x/c) distribution is required to locate 
cp cp 

the chordwlse position of the net span load for constraining the pitching moment 
and is generally selected from a desired two-dimensional section loading. 

The method was applied to the close-coupled wing-canard model of Gloss and 

the results with and without constraints on C are presented in figure 5. 

m 

The variation of the induced-drag efficiency parameter e is a function of 
wing-canard span ratio ^canard/^’ ''^’^8~canard separation h/b, and wing-canard 
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lift ratio. The left side of figure 5 gives estimates of e when is con- 

strained and C is unconstrained. It can be seen that the estimated valu>^.^ 
fn 

of e increase as h/b and Increase and are equal to or greater 

than 1.0. The right side of figure 5 gives estimates of e when both and 

C are constrained to produce a trimmed configuration. It can be seen that the 
m 

extra constraint lowers the values of e; however, e can still be greater than 
1.0 if an upload on the canard is required for trim. When a canard or empnr'nag'' 
download is required for trim, e is equal to or less than 1. In this figure, 
the location of the moment center was completely arbitrary and was chosen simply 
to give uploads and downloads on the canard. 


A detailed study of the e 
fctme.J for one configuration ( 


ffect of moment-center loc..*>tion on e was per- 
h/b « 0.09, **canardA * presented 


in figure 6. It can be seen that e is a maximum at a moment-center location 
of about 10 percent c due to the nearly elliptic span loads present for this 
case. As the moment center is moved away from 10 percent c, the loads required 
on the wing and canard for trim become more nonelllptic and e decreases 
accordingly. 


It should be noted that the wing and canard must be twisted and cambered to 
produce the span loads required to approach the minimum C_ , . The data of 

D, X 

Gloss (ref. 5) were obtained for both flat and cambered wings in the presence of 
a canard. The cambered wings were designed to lift coefficients of 0.35 and 
0.70. These experimental data are compared with the theoretical minimum value 


- S ■ S.o , - 

Qf ijj figure 7. 

theoretical minimum at low 


The uncambered wing alone does not approach the 
because the sharp leading edge does not carry 


any leading-edge thrust. This wing departs drastically from the minimum at 
higher because of the flow separation from the sharp leading edge. The 

downwash and vortex from the canard and strake retard the two-dimensional type 
of separation on the wing and the data show large improvements over the wing 
alone at higher C^, However, the flat wlng-canard-strake combination still 


does not approach the theoretical minimum because of the zero leading-edge 
thrust associated with the sharp leading edge. The cambered wings for the wing- 
canard configuration do approach the theoretical minimum at the design 

because the cambered airfoil carries thrust on the camber line and the leading 
edge is drooped into the local flow direction to reduce the leading-edge flow 
separation. 


Prediction of Induced Drag Due to Control Deflections 

In addition to using the theory to aerodynamlcally design a configuration 
to meet the primary mission requirements, it is also useful in examining the 
effects of deflecting control surfaces and hlgh-llft devices on the Induced 
drag. A modification was made by Paulson and Thomas to the Induced drag mini- 
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mlzatlon program to calculate the Induced drag for specified span loads. The 
Input span loads may be obtained either theoretically or experimentally. An 
example of the variation In span load due to two different types of flaps Is 
shown In figure 8. This analysis was done on the general aviation model shown 
In figure 1 without the fuselage. The span loads were calculated by using the 
Rockwell-Tullnlus method for the plain wing and for the wing with either slotted 
flaps or Fowler flaps deflected. Figure 9 shows the experimental drag polars 
for the three configurations. At C » 1.0, the calculated differences In 

Li 

Induced drag between the plain wing and the wing with slotted flaps or Fowler 
flaps were 0.0010 and 0.0126, respectively. (See table 1.) The corresponding 
differences In the experimental data were 0.0012 and 0.0165, respectively. The 
additional skln-frlctlon drag for the deflected Fowler flap (ref. 6) was esti- 
mated to be 0.0024. When this Is combined with the computed Induced drag, a 
total theoretical Increment In drag of 0.0150 is obtained for the Fowler flap. 
This agrees well with the experimentally measured Increment of 0.0165. 

CONCLUDING REMARKS 


Three applications of theoretical methods for preliminary aerodynamic 
design have been discussed. These methods are used to estimate wing and empen- 
nage geometries and locations to meet performance and stability requirements, to 
estimate span loads for minimum trimmed induced drag, and to analyze the effects 
of control surface deflection on Induced drag. The theories are. In general, 
easy to use, fast, and the agreement with experimental data shows that they give 
accurate results. These methods are being used to design complex multiple 
lifting-surface models for experimental Investigations in the Langley V/STOL 
tunnel . 
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Figure 3,- Comparison between vortex-lattice theory and data. Theory from 
Rockwell-Tullnius vortex-lattice method; data from referebce 5. 
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Figure 4.- Stability margin for model of figure 2. 
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Figure 7.- Effects of canard, strake, and wing camber 
on drag due to lift. 
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